import numpy as np
n = 10;
alpha = 1;
sigma2 = np.array([0.82, 0.91, 0.16, 0.18, 0.28, 0.2, 0.04, 0.07, 0.36, 0.97]);
pmin = 0;
pmax = 1.5;
G = np.array([[13.0, 0.47, 0.84, 0.95, 0.2, 0.05, 0.8, 0.09, 0.48, 0.71], [0.54, 29.0, 0.91, 0.8, 0.87, 0.48, 0.51, 0.95, 0.06, 0.44], [0.26, 0.35, 27.0, 0.12, 0.65, 0.67, 0.39, 0.78, 0.47, 0.28], [0.91, 0.77, 0.72, 11.0, 0.59, 0.84, 0.48, 0.18, 0.38, 0.34], [0.17, 0.04, 0.12, 0.08, 21.0, 0.17, 0.07, 0.08, 0.98, 0.24], [0.66, 0.27, 0.77, 0.78, 0.74, 29.0, 0.48, 0.68, 0.22, 0.07], [0.58, 0.07, 0.8, 0.1, 0.55, 0.78, 18.0, 0.93, 0.86, 0.11], [0.87, 0.16, 0.04, 0.84, 0.48, 0.51, 0.45, 11.0, 0.07, 0.89], [0.97, 0.61, 0.48, 0.18, 0.6, 0.62, 0.64, 0.65, 29.0, 0.51], [0.77, 0.14, 0.9, 0.31, 0.79, 0.31, 0.67, 0.6, 0.85, 23.0]]);